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Abstract We present a strategy for estimation of d-level quantum states and for the simple adaption of 
corresponding measurements. The adaption method is inspired by mutually unbiased measurements, but 
it is also applicable in cases for which no complete set of mutually unbiased bases is known. We present 
results of Monte Carlo simulations, that demonstrate the fidelity gain of the adaptive strategy compared 
to a non-adaptive estimation. 

PACS. 03. 67. -a Quantum information - 03.65.Wj State reconstruction, quantum tomography 



1 Introduction 

From the viewpoint of fundamental quantum mechanics 
as well as from the viewpoint of practical applications 
in quantum information science, determining the state of 
a quantum system P3 H] is an important problem, very 
different from the measuring process in classical physics. 
Perfect reconstruction of a completely unknown quantum 
state requires an infinite amount of measurements for var- 
ious non-commuting observables. One thus also needs an 
infinite number of identical quantum systems. However, 
due to the no-cloning theorem |3j 2] they cannot be copies 
of a single system. As a consequence, initial copies, or 
more exactly, identically prepared systems, elXC el valuable 
resource. From this limited resource only a limited amount 
of information can be extracted. 

This contribution aims at an optimised utilisation of 
this resource for the estimation of a pure c?-level quan- 
tum state by single measurements. The benchmark for 
the quality of such a procedure is the optimal average fi- 
delity of a joint measurement [5-7 on all systems of the 
finite sample. However, the needed physical principles for 
such collective measurements are easily realised. Hence 
we will analyse how well this optimal method can be ap- 
proximated by an approach using measurements on single 
copies and an adaptive choice of measurement parameters. 

For qubits such single-copy measurements have been 
described for pure |MTU] and mixed states [IT]. An adap- 
tive method has been reported in [12] for qubits as well. 
There it was possible to come quite close to the theoret- 
ical limit given in [S]. For higher dimensional states sin- 
gle measurement schemes based on maximum likelihood 
estimation [TJ I13| and mutually unbiased measurements 
|14| have already been used in experiments |15H17| . How- 
ever the application of adaptive methods to higher dimen- 
sions is not straightforward. The methods used for estima- 
tion and adaption cannot be easily calculated in general, 



and the numerical expense to simulate these methods ef- 
ficiently turns out to be quite high. Therefore, alternative 
and simpler methods are needed. Before we explain these 
in detail, let us recall the basic problem. 



2 Basic Concepts 

A finite number of identically prepared versions of a d- 
dimensional quantum system shall be given. We have no 
prior information on the corresponding pure state \ip) and 
want to infer it from the results of single measurements. 
To this end, on one hand we have to reasonably estimate 
a state \ij} est ) from already gathered data (estimation). On 
the other hand we have to decide which measurement pa- 
rameters to use for further measurement steps (adaption). 
We restrict ourselves to projective measurements, which 
can be realised in principle. Hence the main question is for 
the appropriate choice of the measurement basis 
This choice can be influenced by previously accumulated 
measurement data only. Therefore, in the beginning we 
can use an arbitrary basis, e.g. the computational basis. 
But after the first measurement we can adapt the next 
basis according to the previous results, conduct this mea- 
surement, and continue as long as there are copies left. 

In the end, we have acquired a sequence of used mea- 
surement bases and the corresponding results, i.e. which 
vector of each basis was really found in the measurement. 
Based on this measurement data we perform the estima- 
tion procedure, which gives us the estimated state \tp cst ). 
To assess the quality of estimation, we can use the fidelity 
F = | (tp\^p cst ) | 2 . In general, it will of course be different for 
different unknown states |^>). But more importantly, due 
to the probabilistic nature of quantum measurements the 
fidelity even differs for different measurement runs with 
the same state \4>)- While we expect from a good estima- 
tor to work equally well for all unknown states, the depen- 
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dence on the measurement results is inevitable. Therefore, 
to arrive at a valid criterion for characterising the adap- 
tive procedure, we have to average over all possible results. 
In practise we achieve this by averaging over sufficiently 
many Monte Carlo simulations of such measurement runs. 

3 Estimation 

The key part of state estimation is transforming the mea- 
surement data into an estimate of the unknown quantum 
state. A useful method for qubits is maximum likelihood 
estimation [TJ [T3] . But due to the growing number of state 
parameters in higher dimensions sooner or later this be- 
comes highly involved. This is equally well a problem for 
an experimenter, who has to find an estimate for an un- 
known state, as for a theoretician, who wants to assess 
numerically the quality of a measurement scheme. The 
latter has to cope with the additional complication that 
one has to average over a statistically significant number 
of simulated measurement runs. 

Therefore we propose another much more primitive es- 
timation method. For qubits this means calculating the 
average of the Bloch vectors corresponding to the mea- 
sured states and using this average as representation of 
a (mixed) estimated stateQ In order to get a pure esti- 
mated state, this Bloch vector can be normalised. This 
can be viewed as a very rough simplification of maximum 
likelihood estimation, because the likelihood function for 
a single measurement is centred at the measured point 
on the Bloch sphere and the likelihood for more measure- 
ments is the product of such elementary likelihoods. Since 
a Bloch vector is nothing more than a representation of 
the density matrix, this is the same as directly calculating 
the average density matrix and choosing the pure state 
nearest to this average matrix. 

This simple averaging procedure can be generalised to 
higher dimensions d. As in the qubit case a geometric rep- 
resentation of the state can be used. One such generalisa- 
tion of the Bloch vector, that allows this averaging proce- 
dure, is the coherence vector |18l [T5], On the other hand 
the average density matrix p can be calculated directly: 
We simply add all found measurement results {|mfe)} in- 
coherently and arrive at 

v = ^J2\ m kX m k\ (i) 

fc=i 

after v measurements. However, in the subsequent search 
for a pure estimated state, a new complication arises, com- 
pared to the d = 2 case, where it was sufficient to nor- 
malise the Bloch vector of p. 

In order to find the pure state in best agreement with 
the measured data, we choose as estimated state \ip est ) 
the one with maximal overlap with the averaged density 

1 Note that this method gives unbiased information on the 
direction of the Bloch vector only, while its length is not only 
determined by the purity of the unknown state but also by the 
distribution of the measured states. 



operator p. For eigenvalues ej and eigenvectors \ej) of p 
we arrive at 

(r st \7>\r st ) = T t e j \ Cj \ 2 . (2) 

with cj = (e 3 -\tp cst ). When we arrange the eigenvalues in 
descending order eo > e± > . . . > e^-i, we find 

(r st \v\r st )<e , (3) 

where equality holds only if Co = 1 and all other coeffi- 
cients vanish. 

This means that the wanted estimate \ip est ) = |e<j) 
is the eigenvector corresponding to the largest eigenvalue 
eo of p (and the eigenvalue eo itself is the overlap of p 
and \ip cst )). Thus the estimation procedure consists of two 
numerically simple steps: incoherent superposition of the 
measured states and solving the eigenvalue problem of the 
averaged density matrix. For a comparison of this estima- 
tion scheme to other possible estimators applied to non- 
adaptive measurement setups see |Z0] . 

4 Least Bias Adaption 

We have described a method for calculating estimated 
states from measurement data; but we still have to de- 
cide which measurement bases to use. In the qubit case 
it was possible to formulate an expected fidelity with re- 
spect to the used measurement parameters and to find the 
optimal parameters by maximising this expected fidelity 
|12l |2"T] . Again this becomes a tedious numerical task in 
higher dimensions. Instead we use a method that mimics 
mutually unbiased bases (MUBs) [14 . 

Two bases {\bi}} and {|^)} are called unbiased if each 
pair of basis vectors from different bases has the same 
overlap 

\Wi)\* = b ( 4 ) 
In d dimensions there are up to d + 1 bases that are mutu- 
ally unbiased [22 . Such sets of MUBs are optimal for state 
estimation with fixed measurements on single copies in the 
sense of minimising statistical errors [141 [T7] . However, as 
soon as there are more measurements performed than are 
MUBs known, it is possible to find measurement strate- 
gies that improve estimation fidelity compared to repeated 
MUB measurements Our least bias adaption method 
broadens the idea of mutually unbiased measurements for 
use with more measurements. 

The ultimate goal of adaption is finding measurement 
parameters that maximise the information gained in this 
measurement. Mutually unbiased measurements ensure ex- 
actly that, because each sequence of possibly measured 
state vectors is "unbiased", having overlap h with each 
other. However, since the measurements are made succes- 
sively, mutual unbiasedness is sufficient but not necessary 
for this maximised information gain. Only the actually 
measured basis vector has to be remembered from each 
measurement; the directions of the other basis vectors are 
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irrelevant. Therefore the next measurement basis has to 
be "unbiased" with respect to the previous measurement 
results but not to the previous measurement basis vectors 
which were never actually found in a measurement. Us- 
ing this relaxed definition of unbiasedness, more measure- 
ment bases can be constructed. As another consequence 
of this concept of unbiasedness, it depends on the previ- 
ous results which measurements can be called unbiased. 
Therefore, this is a truly adaptive method. Of course, at 
some point Hilbert space is "filled" and no more unbiased 
bases can be found. Then the adaption method must be 
changed to not choosing an unbiased basis, but rather the 
basis with least bias with respect to the measured vectors 
{\m k )}- 

For developing this idea into an adaption algorithm, 
we have to define a function assessing bias in the above 
sense. It should yield unbiased bases as long as this is still 
possible. A good choice for such a function is the entropy- 
like function 



\ 2 \n\(m k \b 3 )\ 2 . 



(5) 



j=0 k=l 



where {|rn/c)} are the v vectors corresponding to already 
found results of the v preceding measurements and 
is the orthogonal basis to be adapted. 

This can be interpreted as the conditional entropy h = 
J2j,kP( m k)p(bj\mk)lrip(bj\m k ), if we set all p(m k ) to be 
equally "true" and thus equally likely [231 [23]. Hence h 
quantifies the information gained when measuring in { \bj) } 
if the "signals" {|mfc)} are already known. 

The basis maximising conditional entropy will be adap- 
ted and chosen for the next measurement. This will be 
done by the equal distribution |(mfc|6j)| = h Vfc, j, if this 
is possible, i.e. there are not too many measurements, and 
we thus find a basis unbiased with respect to the {|wifc)}. 
At some point there will be too many terms in equation 
(|5| to allow this kind of unbiasedness, but the adaption 
algorithm can still choose the basis maximising h for the 
next measurement. 

In order to maximise the conditional entropy |5| cor- 
rectly, we have to ensure that the set is a proper 
basis, i. e. the elements of are not arbitrary vectors, 
but have to fulfil the condition of orthonormality. Since all 
bases can be generated from an arbitrary computational 
basis by a unitary transformation, this is equivalent to 
parameterising a unitary matrix. We did this using the 
Hurwitz construction |25l 126] , producing such a unitary 
matrix from elementary two-dimensional rotations. 



the matrices, we drew their parameters according to the 
relevant Haar measure |25l I2S] . 

Each estimation run uses up to 50 measurements (and 
would therefore require the same amount of identically 
prepared systems) . Since the probabilistic nature of quan- 
tum measurements heavily influences the fidelity of each 
single run, we average over a statistically significant amount 
of such estimation runs (typically 10 4 ). This average fi- 
delity is shown in the following figures with respect to the 
used copies. 



5.1 Qubits 
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5 Monte Carlo Simulations 

To demonstrate the influence of our least bias strategy 
on estimation fidelity in different dimensions, we present 
Monte Carlo simulations of state estimation procedures 
using least bias estimation and compare these to a non- 
adaptive estimation scheme. The latter is realised by choos- 
ing the unitary matrix that rotates the computational ba- 
sis randomly. In order to arrive at an equal distribution of 



number of measurements N 

(b) fidelity difference 

Figure 1. Least biasedness adaption for qubits. The blue tri- 
angles denote least biasedness adaption, using equation Q, 
compared to random measurements (red squares) and the the- 
oretical maximum (grey line). In (a) the absolute estima- 
tion fidelities are shown, whereas (b) depicts the difference 
AF = F — F opt between the simulation fidelities F and the 
optimal fidelity |6|. 
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Figure 2. Least biasedness adaption for different dimensions. The fidelity difference AF = F — F opt to the optimal fidelity ^ is 
shown for four exemplary dimensions d. The blue triangles denote least biasedness adaption compared to random measurements 
(red squares). 



The described method can be used for adaptive estima- 
tion of qubits. The result of the corresponding simulations 
is depicted in figure [l] 

The resulting estimation fidelities are close to the the- 
oretical maximum 

Fopt = w+2 W 

of collective measurements [5J [B]. Furthermore they are 
in good agreement with results reported in |12| . although 
here much simpler methods are used for adaption and es- 
timation. 

5.2 Qudits 

The real test for the presented methods is their applica- 
tion to higher dimensional states. In figure [2] we show the 
estimation fidelities for states of dimension four, six, eight 
and 13. 



In all these examples, there is a gain in estimation 
fidelities due to adaption compared to the non-adaptive 
procedure. On the other hand, in contrast to the qubit 
case, there is a significant gap between adaptive single 
measurements and theoretical maximal fidelity 

p — H±l (7\ 

-Topt — N + d \l ) 

in d dimensions [7]. 

The case d = 6 is of particular interest. Since 6 is not 
a prime power, there is no known way to construct the 
seven MUBs necessary for a Wootters-Fields measurement 
strategy [T3]. In fact there are only three MUBs known 
for d = 6 p7] . But least bias adaption manages to find six 
measurement bases unbiased with respect to the measured 
vectors. Again, this does not mean that the measurement 
bases are mutually unbiased, but that they have an over- 
lap of i only with the basis vectors corresponding to the 
really observed measurement outcomes. Nonetheless the 
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adaption fidelities are the same that would be achieved 
with six fixed MUBs. For the seventh and all following 
measurements the adapted bases are no longer "unbiased" 
but have only least bias with respect to the previously 
measured states. 



6 Conclusions 

We have shown that simple adaptive measurement meth- 
ods can be successfully applied to higher dimensional quan- 
tum states. We presented a method for estimating and 
adapting in higher dimensions which reconciles the con- 
cepts of adaption by information entropy maximisation 
and mutually unbiased bases. 

The simulated estimation fidelities using these meth- 
ods, although smaller than the theoretical maximum, show 
a clear fidelity gain of up to 5% compared to non-adaptive 
methods. This is independent of the dimension and there- 
fore also works for dimensions for which no full set of 
mutually unbiased bases is known. 
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